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Abstract. A new mathematical ansatz is developed for solution of the time-dependent 
Ginzburg- Landau nonlinear partial differential equation describing metastable state relaxa- 
tion in binary (solute + solvent) non-critical solutions with non-conserved scalar order 
parameter in presence of a gravitational field. It has been demonstrated analytically that 
in such systems metastability initiates heterogeneous solute redistribution which results in 
the formation of a non-equilibrium singly-periodic spatial solute structure in the new 
solute-rich phase. The critical radius of nucleation and the induction time in these systems 
are gravity-dependent. It has also been proved that metastable state relaxation in vertical 
columns of supersaturated non-critical binary solutions leads to formation of the solute 
concentration gradient. Analytical expression for this concentration gradient is found and 
analysed. It is concluded that gravity can initiate phase separation (nucleation or spinodal 
decomposition). 


1. Introduction 

Experimental study of the metastable state relaxation in binary solutions is of great 
scientific and industrial interest. In these experiments the usual routine is to investigate 
different stages of the growth process of stable solute nuclei and phenomena associated 
with this growth [1-3]. Metastable state relaxation is not only the growth of nuclei 
which is very rapid but also solution ordering. In this relaxation process very little is 
known about the metastable stage preceding nucleation (mspn). mspn manifests itself 
in solute fluctuations and redistributions which lead to formation of the critical solute 
nuclei. The mspn lifetime is known as induction time t c [1-3]. Recently, experimental 
studies of mspn in vertical columns filled with supersaturated solutions have provided 
interesting data on solute sedimentation [4-8]. In these experiments, initially undersat- 
urated solutions (usually aqueous solutions) have been subjected to a temperature 
quench process which transfers the system to a metastable state. Having carried out 
this transfer it is possible to observe a solute sedimentation process which manifests 
itself in formation of the solute concentration gradients along vertical columns filled 
with the supersaturated solutions. It has been observed that the solute concentration 
gradient appears in vertical columns only at the temperatures corresponding to saturated 
and supersaturated solutions. These experiments dealing with formation of the solute 
concentration gradient in mspn are of great interest since they help to provide informa- 
tion about the nature of mspn by examining solute sedimentation. 
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In this introductory section let us describe qualitatively physics of the solute 
sedimentation phenomenon which takes place in vertical columns filled with supersatur- 
ated binary solutions. Metastable state relaxation is characterized by an induction time 
t c . This time depends on the depth of penetration into metastable region and on the 
solution purity: in the case of little dust and dirt, induction time can be significantly 
increased. In this situation it is possible to investigate different kinetic phenomena in 
the mspn. The only restriction of such investigations is associated with the characteristic 
times of the process under consideration. For example, in our case of the solute 
sedimentation process its characteristic time (sedimentation time) t s has to be: 

(a) less than or equal to the induction time / c ). 

(b) sufficiently long for a noticeable solute concentration gradient to appear. 

At present two attempts to explain this phenomenon have been made [9, 10]. The first 
one [9] proves that no reasonable value for t s can be obtained within the formalism 
of classical non-equilibrium thermodynamics. It has been concluded in this paper that 
the gravitational and buoyancy forces cannot alone be responsible for the solute 
sedimentation process. In paper [10] we have concluded that this process can be 
understood as the redistribution of the solute subcritical nuclei in the presence of a 
gravitational field, mspn is characterized by the situation when new-born solute aggre- 
gates (nuclei) are subcritical, i.e. it is energetically favourable for these aggregates to 
dissolve rather than to grow. In its turn, as it has been demonstrated in [10], the 
probability for these solute subcritical nuclei to appear depends on the density of their 
potential energy hg in a gravitational field, where h is the height of the location of 
the solute aggregate under consideration and g is the free-fall acceleration. Thus, it 
has been derived in [10] that due to the non-linear character of interaction between a 
gravitational field and the field variable describing solute concentration, the probability 
of birth of the subcritical solute nucleus on column bottom is greater than the same 
probability at the top. The only restriction imposed on the theoretical approach 
presented in [10] is that the supersaturated solution has a metastable state which is 
close to the coexistence (binodal) line. In the present paper this restriction will be 
eliminated. 


2. Statement of the problem 

In this paper we will consider an inhomogeneous irreversible process describing the 
relaxation from a non-equilibrium metastable state to an equilibrium stable state in 
the presence of a gravitational field. Equations describing such relaxation processes 
in time are known in the literature [11, 12] as the time-dependent Ginzburg- Landau 
(tdgl) equations. They are formulated in terms of the local order parameter field 
(p(r, /), where r is its spatial location and t is the current time. In our problem, the 
role of this order parameter field can be played by the averaged local solute concentra- 
tion in the new solute-rich phase. At a time just after a temperature quench (transfer 
into a metastable state) there is yet no new solute-rich phase. However, since the system 
under consideration is already in a metastable state there is a stochastic birth process 
of subcritical solute-rich nuclei. These nuclei born due to thermal fluctuations are 
supposed to dissolve in the course of the relaxation process until the moment when 
spontaneous nucleation takes place. At that moment the characteristic size of the 
nucleus (its surface radius, for example) becomes equal to or greater than the critical 
size r c . The time characterizing this moment is the induction time t c mentioned 
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previously. At that time the probability of finding at least one critical solute-rich nucleus 
becomes equal to unity. Thus, induction time characterizes the moment when critical 
and supercritical solute-rich nuclei begin to appear in a solution. These nuclei are 
supposed to grow in the course of the relaxation since it is energetically more favourable. 
Thus, the order parameter <p(r 9 t) corresponding to solute concentration in the new 
solute-rich phase is subjected to small positive changes around zero due to subcritical 
nuclei until the moment of spontaneous nucleation when it begins to grow rapidly. 
This means that the order parameter is a non-conserved quantity with respect to 
metastable state relaxation: 


d t 


d 3 r <p(r, t) 7* 0 
v 


where V is the solution volume. In addition we have to note that for the case of the 
binary solutions order parameter, <p(r , /) is a scalar function. 

Simple theoretical consideration of the situation described above for the spherically 
symmetrical case deals with concept of the minimum work F min (r s ) needed to form a 
solute nucleus of radius r s : 

Emin (O = £> s 2 + E h rl ( 1 ) 

where the first term describes the energy increase associated with the positive energy 
change due to surface formation and the second term is associated with the change 
of bulk energy. It is assumed also that there exists the distribution function W(r s ) of 
solute-rich nuclei with respect to their size r s . Generally speaking, this function depends 
not only on the thermodynamic state of a solution but also on the nuclei flux J . It is 
convenient together with distribution function W{r s ) to introduce the pseudoequili- 
brium distribution function W^rJ for the hypothetical case when there is a detailed 
thermodynamic equilibrium between nuclei dissolution and growth processes (in this 
case J - 0). There is the following expression for function W eq (r s ) [11, 12]: 

We q (r B ) = const -cxp{-E min (r s )/kT} (2) 

where T is the solution temperature. On the basis of these results obtained within the 
classical non-equilibrium thermodynamics it has been concluded by Volmer, Weber 
[13] and Frenkel [14] that the metastable state lifetime / c is proportional to the following 
Boltzmann’s factor: 


1 

,c_ W(r e ) 


W'(r c ) = const • exp{-£ min (r c )//cr}. 


(3) 


It is understandable that the above presented discussions are valid only for the 
homogeneous state. However, taking into account the influence of a gravitational field 
on the metastable state relaxation requires us to consider a heterogeneous case. This 
means that all characteristics of the metastable state relaxation process such as induction 
time t c and critical radius r c become dependent on the density of nuclei potential 
energy hg in a gravitational field. 

A metastable state relaxation equation in terms of the order parameter <p(r, /) was 
introduced by Ginzburg, Landau [15] and Cahn, Hillard [16]. This equation, known 
at present as the tdgl equation, has the following form: 


d<p(r, t) ^ 8Fy(t) 


= -r 


8(p(r , /) 


+ r(r, t) 




rF[<p(r, /)] 


dt 


(4) 



2712 


A F Izmailov and A S Myerson 


where I’>0 is the so-called Landau-Khalatnikov damping coefficient which sets the 
time scale of relaxation process [15] and functional F[<p(r, /)] is the Gibbs free energy 
density for the configuration characterized by the order parameter <p(r , t). The random 
function y (#*, r) represents the Gaussian thermal noise which satisfies the fluctuation- 
dissipation theorem [17, 18]: 

<y(r, f)> = 0 (5a) 

<t)y(r 2 , l 2 )) = ^-^r ^S(r,- r 2 )g(f|- t 2 ) (5b) 

where V r is the volume in the reciprocal space, i.e. in the space of Fourier images. 
Thus, by relations (4) and (5a, 5b) we have determined the Langevin equation. In 
terms of this equation it is possible to distinguish the slow and fast changing degrees 
of freedom: the order parameter <p(r, t) is the slow-changing hydrodynamic degree of 
freedom whereas the random function y(r, t) can be treated as the fast-changing degree 
of freedom which plays the role of a thermal bath. This means that the fast-changing 
degrees of freedom come to almost immediate equilibrium and are described by the 
Gibbs distribution for any given configuration of the slow-changing hydrodynamic 
degree of freedom < p(r , f). The latter, according to equation (4), relaxes toward the 
stable state via the formation of nuclei of the new solute-rich phase. 

Since the Gibbs time-dependent free energy functional F^i t) in equation (4) does 
not depend on relaxation dynamics we adopt, for its density, the standard Ginzburg- 
Landau expression [15]: 

FW(r, f)] = |[V<p(r, t)] 2 +f[<p(x, t )] (6) 

where the coefficient K gives the range over which spatial inhomogeneities persist and 
the functional f[<p( r, /)] is the uniform potential density in the field of which the scalar 
order parameter <p(r, t) evolves. Let us consider the case of non-critical solutions, i.e. 
solutions which do not possess critical points or such solutions where critical points 
are hardly achievable. Partial thermodynamic equilibrium at metastable states of these 
solutions can be characterized at any time instant by the time-dependent Gibbs free 
energy functional F^(/). The simplest form for the potential density f[<p(r 9 /)] of this 
functional in the case of non-critical solutions in a gravitational field can be presented 
as follows: 


r)]=“^ 2 (>, ( 7 ) 

where /x and 77 are the positive coefficients characterizing depth of penetration into 
metastable state, D — VK is the diffusion coefficient and hg is the density of potential 
energy in the gravitational field for the given order parameter <p(r, f) configuration 
(positive direction of the column height is from the column top to its bottom). Now 
one can rewrite equation (4) in the form of the nonlinear inhomogeneous partial 
differentia] equation: 

d<p(r, t) = f) — r Ms p(r, t) + FTi<p 2 (r, t) + y(r, O + (8) 

This equation can be treated as one of the possible modifications of the general tdgl 
equation [15]. The simple analysis of metastability indicates that it can be obtained 
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only for the order parameter <p(r , /) values from the following interval: 

<p s [l-(l -f(/ig)) ,/2 ]<<p(i-, 0<«P* = = !• ( 9 ) 

2 77 I Dfi 

It is noteworthy that gravity can induce spinodal decomposition when hg — T DfjL 2 /4rj. 
In this case, width of the metastable region is negligibly small and any temperature 
quench brings the system directly to the region of unstable states which immediately 
phase separate. For the order parameter <p(r, t) values which are greater than <p s we 
also obtain the solution unstable states. These states cannot exist in nature like the 
uniform states and phase separate through the spinodal decomposition mechanism. 
In addition, result (9) gives us the natural restriction imposed on parameters of the 
tdgl equation describing metastable state relaxation in a gravitational field: 

f(Ag)« 1- (10) 

Thus, with the help of equation (8) one can describe metastable state relaxation 
in a gravitational field for the case of binary non-critical solutions. One of the main 
goals we are striving for in this paper is to determine how the metastable state lifetime 
t c depends on gravity: t c = t c (hg ) in the particular case when solution supersaturation 
is kept unchanged up to the moment of nucleation onset. Exactly this situation has 
been investigated in experimental studies [4-8]. 


3. Solution of the metastable state relaxation problem 


In this paper we are developing a new ansatz for solution of nonlinear inhomogeneous 
partial differential equation (8). For this purpose let us rewrite equation (8) in the 
equivalent form of the following system of two differential equations: 


dipjr, 0 

dl 


= 2 D 


1 A(hg) 


d<p(r,l) i hg i , 

— i y ( r < ') 

dr D 


DV 2 tp(r, I) -2D 


rjhg) J 
\__Mhg) cUpjr, t] 

L r r c (/ig)J dr 


+ I><p(r, l)-r?|f(r, f) = 0 


(11a) 


(nb) 


where r c = r c (hg) is the gravity-dependent critical radius of nucleation which is time- 
independent in the problem of constant supersaturation. The gravity-dependent 
auxiliary function A{hg) will be defined later. 

The mathematical ansatz for solution of system (11a, 116) of differential equations 
can be given in the form of the following scheme: 


(1) We start solution of equation system (11a, 116) from the solution of equation 
(116). Let us try to find the spherically symmetrical solution of this equation in the form: 

<p{r,t) = <p(r,t) = ocx[r(t)] r(f) = |^r-t 


where a and A are the constants to be determined and r s (t) is the time-dependent 
radius of the evolving embryo surface. Substitution of this expression for the spherically 
symmetrical order parameter <p(r 9 t) into equation (life) allows one to rewrite it in the 
form: 


d^rtKO] , , 

dr(() 


A(hg) dx[K0] 
r c (hg) d r(l) 


Tu 

+ -gx[r(t)]~ 


F Tja 
~~D 


* 2 [r(0] = 0. 


(13) 
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As it has been shown by Painleve [19] this ordinary nonlinear homogeneous differential 
equation of the second order can possess solution only under the following conditions: 


2^=5 d 

r c {hg) 



T 7)a 
D 


= 6 . 


(14) 


These conditions allow us to express gravity-dependent critical radius of nucleation 
r c (hg) together with constants d and a through the gravity-dependent function A(hg) 
and parameters (x and rj of the potential energy density f[<p{r , t)]: 


2 I6D 

rAhg) = -A(hg)yj — 

a _6P_25/M\f r c (ftg) r 

“ fi7 4 \T?/L-4(^g)J 


2A(hg) _ /!> 

5 r c (hg) V 6 D 


(15) 


Under conditions (14) equation (13) acquires the following form: 


i 2 *[K/)1 

r\t) 


d ^ (f)] + 6d >[r ( f )]-6^ 2 [r(r)] = 0. 
dr( t ) 


(16) 


Solution of this equation was given by Painleve [19]: 


Af[rU)] = (dC,) 2 e~ 2dr< ' ) P[C, e~ dr<,) + C 2 ; 0, — 1]. 


(17) 


In this expression P[C, e“ dr(,, + C 2 ; 0, -1] is the doubly-periodic Weierstrassian elliptic 
function (dpwef) with invariants g, = g 2 = 0 and g 3 = — 1 (concerning dpwef see Appen- 
dix and [20]). Arbitrary constants C, and C 2 of solution (17) have to be determined 
from the following initial and boundary conditions, respectively: 


lim*[r(/)]| I .-r,(»> = 0 (18a) 


.. d*[r(f)]|„ 

lim — — 

.-o dr s (f) 


rjt ) 


= 0. 


(186) 


The meaning of initial condition (18a) is quite natural: there is no solute-rich phase 
just after the temperature quench process which brings the system into a metastable 
state. In its turn, boundary condition (18 b) implies that solute concentration changes 
sharply on the boundary between solute-rich and solute-poor phases. Thus, conditions 
(18a, h) help us to conclude that: 


C, = z 0 exp 


h'T-fc) JJ 


C, = 0 


(19) 


where r s (0) is the characteristic size of solute unit (molecule or atom) and z 0 ~ 
0.761 647 9998 is the dpwef P[z; 0, -1] zero (for details see in Appendix (A. 16, A.17)): 


P[z () ; 0,-1] = 0. 

Now solution (17) can be rewritten as follows: 

*[r(f)] = (dC,) 2 e ■ 2dr,, 'P[C 1 e- drt, ’;0.-l]. (20) 

where constant C, is given by expression (19). Expression (20) for solution of equation 
(lift) is our final result of the first step in our ansatz for the solution to equation (8). 
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The most interesting conclusion concerning this solution is that it describes appear- 
ance of the gravity-dependent damping spatial periodic structure in the new solute-rich 
phase during the process of metastable state relaxation. It means that this relaxation 
occurs in such a way that solute is redistributed periodically from the very beginning 
of this relaxation process. The time- and gravity-dependent period Ar(f) of such spatial 
heterogeneous structure for solute density can be found from the relation: 


C, expl -d 


4rl,) 4^) 


])■ 


2 w 7 


where oj 2 ~ 0.952 184 9997 is the only real half-period of the dpwef with invariants 
gi = gi = ® and g 3 = — 1 (for details see Appendix (A.1-A.9)). From this relation it is 
straightforward to obtain that: 


Ar(f) = r s (0) + ^ln(^j- 


} 1 raO-r s A (0) 

a e c -'(h g ) ■ 


( 21 ) 


In order to analyse the consequences which follow from this expression for the time- 
and gravity-dependent period of the solute spatial distribution in the new solute rich 
phase it is natural to assume that initially, at the time instant t = 0, solute is in the 
form of elementary solute units with characteristic size r s ( 0). Taking into account this 
assumption together with equation (21) one can conclude that the spatial period of 
the initial solute density distribution in metastable state is gravity-independent: 

Ar(0) = r # (0) + -ln|-pj “i ln(l)r c (0) * 2.290 726 83^(0). (22) 


In this expression we have assumed that r s ( 0)« r c (0) and for values of co 2 and z {) see 
in Appendix relations (A.9) and (A. 17), respectively. This means that initially, just 
after a quench of a solution into a metastable state, all solute elementary units (for 
example, molecules or atoms) are already redistributed in such a way that their spatial 
density is gravity-independent and changes periodically over the entire solution volume. 
This periodic structure of solute distribution is inherent to metastable states of non- 
critical solutions: it is initiated by transition into metastable state and is preserved 
during the following metastable state relaxation process. For the stationary case 
corresponding to formation of the critical solute nuclei the critical period Ar(/ C ) of 
their spatial distribution becomes gravity-dependent: 


* \ , mi5 , a , 1 r c (hg)-rJ0) 

M/ c ) = r s (0) + 2 ln( 2 )r c ( 0 )+- ~ lf - - ■ — . 

A r c {hg) 


(23) 


Analysis of this expression will be given later on. 


(2) Now let us take into account agreement between two differential equations (11a) 
and (lib). Thus, substituting spherically symmetrical expression for the order 
parameter <p(r, f) found from equation (lib) into equation (11a) one can obtain the 
following evolution equation for the nucleus surface radius r s (/): 


dr s (t) 

dr 



Mhg) 

.Mbg) 


1 

r 


hg 

aD 


K-'[r{t)] 



A - 


+ P(r , t) 


(24) 


p(r, l) = - 


rjfa) 

. r s (0 . 


A -1 

(aK[r(t)]) ' y(r, t) 


K[r(l )] 


dx[r(t)] \ 

MD ) 


where 
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In order to investigate evolution of the spherically symmetrical boundary of solute 
nucleus it is necessary to rewrite equation (24) for r=r f (l): 


d rjj) 

dt 



Mhg) 

.rjhg) MO 




+ p[M 01 


(25) 


where 


p[MO] 


rc(fcg) 

. MO . 


A - l 

(«K[MO]) 


'y[M01 


K[MO]= K0(0]|r=r, ( n- 

Thus, we have obtained the Langevin equation for evolution of the nucleus surface 
radius MO- Random force p[M01 of this evolution process has the following corre- 
lators: 


(p[MO]> = 0 


<p[M'.)]p[MO)]> = D[r,(r,)]«(/,-i 2 ) 


where 


D[MO] = T 


r JM) 

. MO . 


2( A 


(«K[MO]) 2 


(26a) 

(26b) 


(27) 


and averaging is over all possible realizations of the random force p[M 01- 

Let us now define auxiliary function A(hg) from the condition that at stationary 
state there exists the following equality: 


lim 

r r t 



(28) 


From this condition one can obtain that: 


A{hg) = 1 + 


hgrJJ ig) 
2aD 2 


K-'[r c (bg)] 


A ( 0)=1 


(29) 


where 


K[r c (hg) ] = lim X[r s (0] = — d(dC,) 2 (2F[C, ; 0, -l] + C,P <n [C, ; 0, -1])+ O(e) (30) 

t-*t y 


P ,U [C 1 ;0,-1] = 


dP[z; 0, — 1] 
dz 


z = C\ 


Now taking into account that function /^[r s (/)] is the decreasing one with respect to 
r s (r) it is possible to conclude that: 

(a) when r s (/)< r c (hg), then the averaged nucleus surface radius <r s (r)) is decreasing 
in the course of the metastable state relaxation d<r s (/))/df < 0; 

(b) when r s (t) = r c {hg), then the averaged nucleus surface radius <r s (/)> is kept 
unchanged d(r s (f))/df = 0 until eventually spontaneous nucleation takes place due to 
appropriate random thermal fluctuation; 

(c) when r s (t)> r c (hg), then the averaged nucleus surface radius (r s (t)) is increasing 
in the course of the metastable state relaxation d(r s ( t))/dt > 0. 

The last conclusion concerning growth of the averaged nucleus surface radius r s (t) 
is due to the nucleation process only and has nothing in common with the coalescence 
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process which we do not consider in this paper. Thus, the first term on the right of 
equation (25) determines the systematic force driving the growth (decay) of solute 
nucleus of the new solute-rich phase. This force tends to conserve the system uniformity 
and by itself cannot induce transition into a new phase. Such a transition can be 
induced only by the second term on the right of equation (25), i.e. by the random 
force p[r s (7)] which describes thermal fluctuations. 


4. Fokker-Planck formalism and gravity-dependent induction time 


Let us now introduce a statistical description of the system which is governed by 
equations (25)-(27). A metastable state produced when the equilibrium line of the 
first-order phase transition is crossed with finite velocity is the state of non-complete 
equilibrium. In such a state, distribution of the small-scale degrees of freedom (r« 
r c (hg)), such as thermal fluctuations, corresponds to the local equilibrium. The large- 
scale degrees of freedom (r> r c (hg)), such as the scalar order parameter field <p(r, t ), 
vary in the course of time in correspondence with equation (11a). The critical nucleus 
radius given by expressions (15) and (29) depends on the depth of penetration into 
the region of metastability of the initial phase and on gravity. It has been demonstrated 
in expression (29) that the last dependence is highly nonlinear. 

Metastable state relaxation will be described in this section as relaxation of the 
size-distribution function W[r s (f)] for solute nuclei of the new solute-rich phase: 

W[r s (/)] = [ dr s (OG[r s (/),r s (/')]W[r s (/')] (31) 

J r s (0) 

where G[r s (/), r s (/')] is the transition probability. Following the theory of homogeneous 
random processes [21] we introduce this probability in the form: 

G[ r s ( r, ), r s ( t 2 )] = (8[r s (ti)- r 8 ( / 2 )]>- 

It is well-known [21,22] that when random force in Langevin equation is Gaussian 
distributed it is possible to construct the stochastically equivalent Fokker-Planck 
formalism. Within this formalism one is supposed to deal with the linear partial 
differentia] equation with respect to the size-distribution function W[r s (t)]\ 

dW[r s (t)] _ dJ[r s (t )] (32) 

St dr s {t) 

where J[r s (/)] is the density of flux in the size-space. For this density there is the 
following expression in terms of Langevin equation (25) characteristics: 

i[r,(0] = -B[^(0]VV[r s (f)]-D[r s (f)]^L^ (33) 


where 


B[r s (/)] 



1 


A(hg) ~ [ hg K 
»'c(^g)J a D 



1 


Functions B[r 5 (0] and D[r s (/)] in expression (33) have to be differentiable real 
functions with the only restriction that D[r s (/)]>0. The first term on the right of 
equation (33) describes systematic growth of solute nuclei of the new solute-rich phase 
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and is usually called the ‘transport term’ whereas the second one is responsible for 
diffusive growth of these nuclei and is called the ‘diffusion term’. From the point of 
view of the Fokker- Planck formalism function D[r s (f)] plays the role of the generalized 
time-dependent diffusion coefficient. It is important to note that this coefficient goes 
to zero like rf 1 A> (f) (A > 1) in the case when the surface radius r s (f) of solute nucleus 
tends to infinity. It can be also noted that for the particular case described here diffusion 
coefficient is related to the probability that a solute elementary unit (molecule or atom) 
joins the solute-rich nucleus of size r s (f) per unit time. 

The first step in analysis of equation (33) is to find equilibrium distribution function 
W / eq [r s (/)] corresponding to zero flux density J[r s (t)] = 0. This distribution describes 
the partial thermodynamic equilibrium when nucleus growth and decay processes 
compensate each other. Thus, it is straightforward to find the following expression for 
the equilibrium distribution function W eq [r s (0]‘ 

VV'eqt '•»(')]= VV eq [r s (0)] expj^— J df "5(7)]- (34) 


This function gives the equilibrium distribution of solute-rich nuclei with respect to 
their size. 

In addition to equilibrium solution W eq [r s (t)] one can find the steady state solution 
WstWO] of equation (33) which corresponds to the constant stationary flux density 
J[r s {t)] = constant: 


Wstl> s (0]= W eq [r s (f)] 1-J st [ dr - * 

L J,joi D(r)W e Jr) J 


where 


(35) 


J " Lo> dr D(r)VV eq (r) : 


1 


We q [r s (0)] 


r* , i rr b (oi 

dr— ——exp dr — — — . 

J rjo) f7(r) LJr,(o> D(r )_ 


(36) 


The size-independent flux density J sl corresponds to the constant flux density in the 
stationary nucleation regime. It means that magnitude of this density is the nucleation 
velocity, i.e. the average change of the nucleus characteristic size per second. Thus, 
the reverse quantity can be associated with the metastable state lifetime / c : 

t c = J sr (\m)-' (37) 


where 1 m is the unit length. 

The second step in analysis of equation (33) is to find analytical expression for 
function W eq [r 8 (0]* In order to carry out this step let us find the interval [K lu K c ] 
within which derivative K[r s {t)] (for definition see expression (25)) varies when its 
variable r s (/) changes in the interval [r s (0), r c (hg)] (for details concerning derivative 
of the dpwkf see in appendix expressions (A.12 and A. 14)): 

Ko= hm K[r s (r)] = -(dz 0 ) 3 (38a) 

r,< /)-/•„((» 


K c = lim 


KK(f)] = -(dz 0 ) 3 + O 



(386) 


The last equality for is valid only until r c (6g)»r s (0) which is usually the case in 
nucleation problems. From these expressions for K 0 and K c one can conclude that 
derivative K[r s (t)] is almost everywhere constant (/C[r s (t)] = -(dz 0 ) 3 ) when its variable 
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r s (f) changes from r s (0) to r c . Taking into account this result it is possible to obtain 
from relations (15) and (29) the following expression for the gravity-dependent critical 
radius of nucleation r c {hg ): 

r c (hg) = r c m\+f(hg)]-' (39) 

where 


f(hg) = i{hg) =0.61$ 714 749 • {(h g) < 1. 
10zi) 


It is obvious from this result that the critical radius of nucleation slightly diminishes 
in the direction from the column top to its bottom. This decrease becomes more 
considerable for the metastable states close to spinodal line when f(bg) tends to unity 
(see expression (9)). 

Now let us return to the equilibrium size-dependent function We^r^f)]. Taking 
into account results (38a, 38b) it is possible to obtain that there is the following 
expression for the equilibrium distribution of the subcritical (r„(/)=s rjhg)) solute-rich 
nuclei with respect to their size: 


WUr s (r)]= W cq [r s (0)]exp 


a(dz 0 ) 3 


L r r x c -'(hg) 

2Da(dz u 


[ dr 2Da(dz„) 3 r A_2 - 
J r s < 0) L L 

a(d z 0 ) 3 / ar , (J J 

P — I 2Da(dz„) I 


Mhg) A .1 _^ r A ll 
r c (hg) \ D J 

~~ R x '(f, hg)-A(hg)R x (t, hg ) 

A 1 




| bgr c (/ig ) 
D 


R x (t,hg) 


(40) 


where R(t, hg) = r s (t)/r c (/ig) and function A(hg) is given by expression (15). From 
this expression it follows that for the equilibrium state under conditions r 5 (0)« 5r c (0)/2, 
r (/) « r c (hg) and A = 3, it is possible to restore phenomenological expression (1) for 
the minimum work E min (r s ) needed to form subcritical or critical solute nucleus of the 
new solute-rich phase with characteristic size r s : 

E mi „(r s ) = E s r 2 +E„r 3 (41) 


where 


E s = 


E s (/ig) = 


Po 2 (dz u ) 6 

rrl(hg) 


E h = E h {hg) — - 


2Da 2 (dz„) 6 

3Fr 3 (/ig) 


bga(dz„) 3 
3T Drl(hg) 


It is noteworthy that within the developed above theoretical approach we have obtained 
the correct sign for coefficient E h without special assumptions. Usually the first term 
of expression (38) is associated with the surface energy: E s = 4-n-er, where a is the 
surface tension coefficient. Thus, within our approach one can express this coefficient 
through the tdgl model parameters: 


kTa = kTcr( hg) = 


3Pa 2 (dz ( ,) 6 

4v\Vrl(hg) 


75 

576 7rA 



[1 +fihg)f. 


(42) 
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What concerns height h dependence of the surface tension coefficient <x(/ig), it can 
be concluded, that for the fixed gravity acceleration g it grows in the direction from 
the column top to bottom. 

Analysis of result (40) gives that the function W~ q '[r s ( f)] has a sharp maximum at 
r s (f) — r c {hg) which reflects the existence of a nucleation barrier. The expression for 
this maximum has the form: 

VV £q ‘[r c (%)] = e w (43) 

where 


Li(hg) = 


2D« 2 (dz„) 6 


AT 


A — 1 


A(hg) 


hgq{dz a ) 3 r c (hg) 

AT D 


2Pa 2 (dz 0 ) fe ^4 1 

AI'(A-I) _ 3A(A — 1 ) -? 7 2 ~~ r c (0) 


In order to obtain this result we have used expressions (29) for A(hg) and (38 b) for 
K[r c (hg)]. It is important to note that function L A (/ig) = L x does not depend on gravity. 

The third step in analysis of equation (33) is to arrange adequate estimations for 
the induction time t c . In order to carry out these estimations let us calculate the integral 
in expression (36): 


where 


j:.=- 


— - — f 

V^ q [r s (0)]J, 


dr — — 
( 0 ) D{r) 


^Alr) 


(44) 


A(r) 



Mhg) \ 

r,{hg)J 


K(r') + 


hg 

2D 2 a_ 


(r') A 'K(r') 


e 


2D Tj 
fff 3 


It can be done with the help of the steepest-descent method [23] which can give correct 
asymptotic for integral (44) only in the following case: (a) e = rj/3 » 1; (b) functions 
D ] (r) together with A{r) are the real and sufficiently smooth functions of r. The 
classical result for such integrals is due to Laplace who has argued that their main 
contribution comes from the neighbourhood of the global maximum of function A(r). 
This maximum takes place at r = r c (/ig) and assuming now that A(r) is twice differenti- 
able function it is possible to obtain that: 

where 


277 


i) 


M ■ '[ r c ( hg) ] ) D[ r c ( hg ) ] lV eq [ r c ( hg )] 


(45) 


A i2) [r c ( hg)] = 


d 2 A[r s (f)] 

dr 2 (f) I r,( r ) = r c { h% ) 

«’ „ 2 r t , ^ hga 2 dK[r s (t)] 

r c (hg) 2D dr s (f) 


(46) 


r s { t ) = r L ( hg ) 


Dlr c (hg)] = ^K- 2 Mhg)]. 


(47) 



Metastable state relaxation and gravity 


2721 


As it has already been demonstrated there is the following approximate equality: 
K[r c {hg)]^ -(dz 0 ) 3 r s (0)^r s (t)^r c (hg) 

which is correct if, and only if, condition (38 b) is satisfied. Thus, substituting 
expressions (46-48) into equation (45) one can obtain an expression for the gravity- 
dependent induction time t c in the form: 

*c=*c(fcs) = U0)[l +f{hg)Y 1 (48) 


where 


* c 




^ eq [r c (0)] 


*(lm). 


In order to obtain expression for stationary fiux density J sl in form (48) we have used 
expressions (15) and (39) for parameters a, e and function r c (hg), respectively. Relation 
(48) for induction time demonstrates that this time is gravity-dependent. For the case 
of the fixed gravity acceleration it is straightforward to conclude that induction time 
diminishes in the direction from the column top to bottom. Taking into account that 
maximum value of function f(hg) is equal to 3 /IOzq- 0.678 714 749 ■ {((hg) = 1) = 
0.678 714 749 (see expression (39)) it is possible to conclude from result (48) that the 
minimum achievable induction time in nucleation experiments in vertical columns is: 


^min = 0.595 693 8191 • r c (0). 

The following decrease of induction time due to gravity is impossible since the following 
increase of column height at fixed gravity acceleration leads to phase separation on 
the bottom level of vertical column through spinodal decomposition mechanism (see 
expression (9) for boundaries of the metastability region ). Taking into account relation 
(9) for <p s one can rewrite expression (48) for induction time / c (0) in the form: 


f c (0)=-z^ s 


te)' 


1 


WUr c (0)] 


■(lm). 


(49) 


Now it is apparent that induction time t c (hg) depends not only on gravity but on the 
depth of penetration into metastable region. To explain that one can distinguish two 
limiting cases: 

(a) <p s ~ 0. In this case for the fixed value of hg time t c ( hg)~ 0 what corresponds 
to the infinitely narrow metastable region (binodal line separates stable and unstable 
states), i.e. in this case temperature quench brings system almost on the boundary with 
unstable states which immediately phase separate. 

(b) <p s »1. In this case for the fixed value of hg time t c (hg) may be very large 
corresponding to the very narrow unstable region, i.e. in this case the temperature 
quench always brings a system to a metastable state which can be thought to be very 
close to the binodal line (to the region of stable states). 

All these derived conclusions concerning the induction time t c (hg) are consistent 
with the known experimental facts [1-9]. 

By result (41) we have also obtained restrictions under which the well-known 
assumptions (2-3) are correct. These assumptions state that there is the inverse propor- 
tionality between the induction time f c ( 0) and the equilibrium size-distribution function 
W eq [r c (0)] for critical solute nuclei. This conclusion, / c (0) ~ W“ q [r c (0)], was first 
obtained by Zeldovich [24] on the basis of classical non-equilibrium thermodynamics. 
In addition it has been proved that the equilibrium size-distribution function 
We q [r c (/rg)] is gravity-independent: W K J_r c (hg)]= W eq [r c (0)]. 
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5. Conclusions 

All results obtained in this paper are due to a new mathematical ansatz for solution 
of nonlinear partial differential equation (8). This equation describes relaxation in a 
gravitational field of the non-conserved scalar order parameter <p(r, /) which is associ- 
ated in our paper with solute concentration in the new solute-rich phase. The ansatz 
is constructed in such a way that equation (8) has been split into two interconnected 
equations (11a) and (lib). The doubly-periodic solution of second equation (lib) is 
found and analysed. This solution corresponds to the new topology in solute redistribu- 
tion within the new solute-rich phase. According to our result (20) this redistribution 
corresponds to appearance of the damping spatial periodic structure for solute density 
distribution. The developing periodic structure exists during the metastable state 
relaxation process until coalescence processes begin to play an essential role. The new 
heterogeneous topological structure of the solute density distribution is inherent only 
to the solution metastable state: it appears immediately after solution transformation 
into metastable state. Since only one period of the found dpwef is real (see in Appendix 
expression (A. 9)) the solute heterogeneous structure described is in reality a singly- 
periodic structure. The main characteristic of such a spatial singly-periodic structure 
is its period. We have found that this period Ar(f) is gravity- and time-dependent (for 
details see expression (23)). 

It is obvious from solution (20) that the metastable state relaxation is also 
heterogeneous with respect to the density of potential energy in the gravitational field. 
This fact proves the existence of the solute sedimentation process which takes place 
in the course of the metastable state relaxation process in a gravitational field up to 
the moment of spontaneous nucleation onset. Solute sedimentation manifests itself in 
formation of the solute concentration gradient along the column height. This gradient 
is directed downward and is associated with contribution to solute density due to 
solute subcritical nuclei. In order to estimate the solute concentration gradient 
A H <p[r 5 (0] = a{Affo(0]|/,=H-*[r s (')]|/,..o}, in the vertical column of height H let us 
assume that for any stage of evolution for the surface radius r s (t) of solute nucleus 
there exists the following inequality: 



where the surface radius r s (f) changes in the interval [r s (0), r c (hg)]. Underthis condition 
it is possible to expand solution (20) into series around z 0 (see expression (A. 18) in 
Appendix) and to obtain A H <p[r s (/)] in the form: 

A„*[r s (/)] = f^: [r s A (0 - r A (0)][(l +/(//g)) A - - 1]. (51) 

AT C (UJ 

It is apparent from this expression that the solute concentration gradient is an increasing 
function of time and gravity. Analysis of expression (51) for t = t c (Hg) gives that the 
maximum gradient & H <p[r c (Hg)] of solute concentration <p[r s (f)] is the non-trivial 
function of the potential energy density Hg in a gravitational field. It is easy to 
demonstrate that for the case when A = 3 and r s (0)/r c (0)« 1 this maximum achievable 
gradient of solute concentration has local maximum with respect to Hg ~ ( Hg ) max 
5zo(V3 - 1 )YDii 2 /6r) = 0.269 645 9773 • Y Dp 2 fry. 

A„<p[r c (Hg) ma J] s ^|^0-045 368 1232-^. 


(52) 
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Further increase of the potential energy density Hg leads to the decrease of the 
achievable solute concentration gradient since gravity-dependent critical radius of 
nucleation r c ( Hg) and induction time t c (hg ) become sufficiently small for the nucleation 
process to start. Thus, the rapid nucleation prevents formation of the considerable 
solute concentration gradient along the column height. 

It is also important to note that when column height H is sufficiently large: 

rr SzlDVfi 2 
ti — 

6 gV 

it is possible to initiate immediate spontaneous nucleation by gravity. It means that 
gravity-dependent critical radius of nucleation r c (Hg) corresponding to column bottom 
(h-H) becomes equal to the initial characteristic size r s (0) of solute units (molecules 
or atoms). In this case nucleation on the bottom level of vertical columns starts 
immediately after solution transfer into a metastable state. 


r c(0) 

rAO) 


-1 


(53) 


Appendix: Doubly-periodic Weierstrassian elliptic functions 


Doubly-periodic Weierstrassian Elliptic Function (dpwef) P(z;g 2 ,g 3 ) is an even 
function of order two which is usually defined by the following double series [15]: 


P(z; g 2 , 8s) = *>')= <*>') 


= ~2 + Z Z (1 “ 1 

Z m n 




L(z 


J 1 

-O ) 2 o 2 

d L mn / “fflflj 



(A.l) 


where 


il mn = mluj + n 2a/. 

There is also the integral formula for the dpwef P(z; g 2 ,g:0 : 

Z= I d>'(4y 3 -g2j'-g3)~ 1/2 ‘ ( A - 2 ) 

In expression (A.l) w and co' are the dpwef primitive half-periods; m and n are the 
integers. Summation in expression (A.l) ranges over all integers m and n except m =0 
and n = 0 simultaneously. Double periodicity means that there are the following 
expressions: 

P(z + 2o)\(o, co') = P(z\cd, a>') (A. 3) 

P(z + 2(o'\<o, co') — P(z\o>, u) ( ). (A. 4) 

The function P'(z|<u, co') = dP{z\w, (o f )/dz is the odd elliptic function of order three 
with half-periods (o a (a — 1,2, 3): 

CD] — W C0 2 — ~<X> ~ CO (Dy = co' . (A. 5) 

It can be concluded that P’(a) a \<o, (o') = P(— <o (X \<o, (o') — 0 since P'(z\a), co f ) has half- 
periods (o a and P\-(o a \w, w ') = -P'{a> a |w, co). Thus, it follows that z = a>„ (o = 1,2, 3) 
is an irreducible set of zeros for P'(z\<o, co'). It is customary to put that: 

e a = P{(o (Y \(o, co') (a = 1, 2, 3). 


(A. 6) 
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The dpwef invariants g, , g 2 and g 3 can be defined in terms of constants e„ (a = 1, 2, 3) 
in the following way: 

g, = e, + e 2 + e, = 0 

g 2 = 2(e\+e\+el) (A.7) 

g 3 = 4e,e 2 e,. 

Discriminant A of these relations between and g, t (a = 1, 2,3) is given in the form: 
A= 16(e,-e 2 ) 2 (e, -e 3 ) 2 (<? 2 -<? 3 ) 2 = g 2 -27g 3 . 

In our particular case corresponding to g, = g 2 = 0 and g 3 = -l one can obtain that: 

e x =4“ 1/3 e' 7 * 73 e 2 = -4 1/3 e 3 = 4 in e _i7r/ \ (A.8) 


Now taking into account relations (A. 2) and (A. 6) it is possible to write the following 
expression for the real half-period cj 2 : 


(j)-> — VP 


dv(4y 2 + 1) _1/2 = 


5 7T 

3372^571 : 


- 0.952 184 9997 


(A. 9) 


where VP is the Cauchy principal value. For non-zero invariants g 2 and g 3 there exist 
the other representations: 

'g: = 6011(1- 5 m ,„)( 1 - 5„,,)n J, g 3 = 140 1 1 ( 1 - S m0 )( 1 - s„. 0 )n 2 . ( A. 1 0) 

m n m n 


From the definitions presented above one can derive the homogeneity relations for 
arbitrary s ^ 0: 


P(sz|so>, so/) = s 2 P(z |a>, ad) 

P(sz\ s~ 4 g 2 , s _6 g 3 ) = s~ 2 P(z ; g 2 , g 3 ) 


(A.l 1) 


and 


(A. 12) 


P'(sz|so>, sad) = s 2 P'(z |a>, ad) 

P'(sz; s“ 4 g 2 , s" 6 g 3 ) = ^ 3 P'(z; g 2 , g 3 ). 

Taking into account relation (A.l 1 ) and putting s= i one can obtain for the case of 
negative g 3 the following result: 

P(z; g 2 , g.0 = &2 > -g.O* (A. 13) 


Thus, it is always possible to convert the case corresponding to negative g 3 into the 
case corresponding to positive g 3 . 

Zero z ( , of the dpwef is defined through the relation: 

g2,g3) = w') = 0. (A. 14) 


At this point there also exists the following useful relation: 

P’(z <>; g2,gs) = »VgI- 


(A. 15) 


Now taking into account integral formula (A. 2) for the dpwef P(z;g 2 ,g 3 ) we 
apparently find z 0 in the form which corresponds to g 2 = 0 and g 3 = - 1 . We obtain that: 




• A 

dy(4y' 

Jo 


g2.V-gj)'' 


(A. 16) 


dy(4/+l) 



==0.761 747 9998. 


(A. 17) 
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Now it is straightforward to obtain the following expansion of the dpwef P[z\ g 2 , g 3 ] 
around z () : 

P[z-,g2,g 3 ] = P"W,g2,g>Ml~3C 2 u 4 +- • •]-[5C 2 +14C,m 2 +- • •] ( A.18) 

where C 2 = g 2 / 20 = 0 and C,/28=-l/28. 
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